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In this paper we consider the collision of spinning holes using first order perturbation theory of 
black holes (Teukolsky formalism). With these results (along with ones, published in the past) one 
can predict the properties of the gravitational waves radiated from the late stage inspiral of two 
spinning, equal mass black holes. Also we note that the energy radiated by the head-on collision of 
two spinning holes with spins (that are equal and opposite) ahgned along the common axis is more 
than the case in which the spins are perpendicular to the axis of the collision. 
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■ I. INTRODUCTION 

o ■ 

' There is considerable current interest in studying the coUision of two black holes, since these events could be 
primary sources of gravitational waves for interferometric gravitational wave detectors currently under construction. 
(~| ' The mathematics describing such events consists of the nonlinear partial differential equations of general relativity, 
Einstein's theory of gravitation. These equations are very difficult to solve especially in the case of two merging holes, 
or of the highly distorted final hole that is formed by the merger. Limited (in resolution and evolution time) full 
CO numerical simulations of grazing coUisons of black holes are being done now , but long term stable evolutions are 
still distant in the future. It is therefore of interest to have at hand approximate results which in certain regimes could 
be used to test and possibly even complement the full numerical evolutions. Among such approximation methods is 
the "close-limit" approximation in which spacetime of a black hole collision is represented as a single distorted black 
hole and evolutions are done with simple linear perturbative equations. In the past, this method has been applied to 
the case of a head-on collision of two boosted holes and the case of slow inspiral 1^ with considerable success. In 
T-H ' both these cases, the black holes considered had no spin. In this work, we will consider the collision of spinning holes. 
\ Since we are attempting a linear, first order perturbative calculation, it is sufficient to study the coalescence of two 
spinning holes with zero initial linear momentum. This is because, to obtain the results for the inspiral of spinning 
holes, one can simply do a linear superposition of these results with those from our past work ||3|]. However, the 
O i first order perturbative treatment imposes several restrictions on the physical scenarios we can treat. In addition to 
'Y* '' the close-slow limitations, we can only easily analyze cases in which the spins of the two black holes are equal and 
$_H . opposite. This happens because, to lowest order, the (perturbative) extrinsic curvature of a black hole with spin S 
and a conformal distance L from the origin is of order SL. Superposing another hole with parallel and equal spin and 
^ , at —L will yield a situation with zero perturbation in extrinsic curvature, which is not a case of interest. 

So in the following, we shall present two cases that are exhaustive. Case /, when the two spins are equal and 
opposite and aligned along the common axis of the two holes. Specifically, we place the two holes at a; = ±L/2 
I and have their spins parallel and anti-parallel to the x-axis. Case II, when the two spins are equal and opposite 
■ " ' and perpendicular to the common axis of the two holes. Here, we place the two holes at x = ±L/2 and have their 
spins parallel and anti-parallel to the y-axis. These cases have no net angular momentum. So we expect Zerilli and 
Teukolsky formalisms to agree exactly (Kerr parameter a = 0). We, therefore report calculations and results from the 
Teukolsky formalism only. Part of this work, in the context of the Zerilli formalism can be found in this reference [Q . 
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II. INITIAL DATA 



To evolve a spacetime in general relativity, one needs to provide initial data, a 3-geometry gab and an extrinsic 
curvature Kab, that solve Einstein's equations on some starting hypersurface (i.e., at some starting time). These 
initial value equations have the form, 

^''{Kab - gabK) = (1) 

- KabK"^ + ^0 (2) 

where gab is the spatial metric, Kab is the extrinsic curvature and ■^i? is the scalar curvature of the three metric. If we 
propose a 3-metric that is conformally flat gab = 4''^gab, with gab a fiat metric, and 0^ the conformal factor, and we 
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use a decomposition of the extrinsic curvature Kab = (j} ^Kab, and assume maximal slicing K'^ = 0, the constraints 
become, 

V''Kab = (3) 

V^(ly = -\(t>-'kabK''\ (4) 
o 

where V is a flat-space covariant derivative. 



To solve the momentum constraint, we start with a solution that represents a single hole with spin S |10|, 



- ^ [ecad5 + eefcrf^ Vn,] . (5) 



In this expression for the conformally related extrinsic curvature at some point x"", the quantity Ub is a unit vector, 
in the "base" flat space with metric gab, directed from a point representing the location of the hole to the point x"' . 
The symbol r represents the distance, in the flat base space, from the point of the hole to x°' . It is straightforward to 
show that the solution of the Hamiltonian constraint corresponding to eq. (5) corresponds to a spacetime with ADM 
angular momentum 5". 

The next step is to modify this to represent holes centered at a; = ±L/2 in the conformally flat metric. Since the 
momentum constraint is linear, we can simply add two expressions of the above form, 

KTb° = KT {x^x- L/2, or Sy ^ S) + K^^ {x x + L/2, or Sy = ~S) . (6) 

We will choose in further expressions to use a polar coordinate system in the flat space determined by gab centered in 
the mid-point separating the two holes and label the polar coordinates as {R, 9, 0). So i? will be the distance in the 
flat space from the midpoint between the holes. 

To solve the Hamiltonian constraint |4[ we introduce an approximation, (the slow approximation) which we will 
show is enough for our purposes. In fact, in this approximation the solution for the conformal factor turns out to be 
the familiar Misner B solution if one chooses the topology of the slice to have a single asymptotically flat region, or 
the Brill-Lindquist H solution if there are three asymptotically flat regions. 

A. The slow approximation 

We assume that the black holes are initially close, and that the initial spins S are small. We denote by n"*" and 
n~ the normal vectors corresponding, respectively, to the one hole solutions at a; = +L/2 and at a; = —L/2, and we 
recall R to be the distance to a field point, in the flat conformal space, from the point midway between the holes. 
For large R, the normal vectors n+ and almost cancel. More specifically n+ = —n~ + 0{L/R). A consequence of 
this is that the total initial K°'^ is first order in L/ R. It scales linearly with the spin S as well. Thus the source term 
in the Hamiltonian constraint is quadratic in S. If we choose to find a solution to the conformal factor to first order 
in S (which should give us a good approximation in the case of slowly spinning holes), we can ignore this quadratic 
source term. So now, the Hamiltonian constraint looks like the one for zero momentum, which is simply the Laplace 
equation. A well known solution to this, is the Misner solution ||]. This solution, is characterized by a parameter /zq 
which describes the separation of the two throats. We can relate this parameter to the conformal distance L in the 
following way 0, 

zIji ^ — ' smh n/in 

n— 1 ^ 

We must now map the coordinates of the initial value solution to the coordinates for the Schwarzschild (Kerr with 
a = 0) background. To do this, we interpret the R as the isotropic radial coordinate of a Schwarzschild spacetime, 
and we relate it to the usual Schwarzschild radial coordinate r by i? = {y/r -\- \Jr — 2Af )^/4. From this we arrive at 
the following expressions for the components of the metric and extrinsic curvature. 

B. Initial data for the Teukolsky function — Case I 

Following the construction as outlined in the last section, we get the following expressions for the initial spatial 
metric and extrinsic curvature: 
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Kab 



3SL 





-2 sin 9 sin 2yj 



— 2 sin 9 sin 2ip 
^l-2A//r 

r cos sin 2(/7 



-yi-2M/r 

— 2 sin sin 2^^ cos^ 02 • nf 2/1 2 • 2 \ 

, ^ r sm 6^ COS COS 03 — sm 09 



— 2 sin ^? sin 2f? cos ^ 

r sin 0(cos^ 6 cos^ ip — sir? (p) 
—r sin^ 9 cos sin 2(p 



3 (r-2Af)(-l + 2,/l-2M 
4i(cos0±l)2L5 ^ 

4 r(r(l + ^l-2M)_M) 



(8) 



The metric has a form, identical to the Misner solution ||6|. Now, we use the methodology and expressions in ||3|, to 
find the initial data for the Teukolsky function, Vl/ = p~'^ip4, where p = — l/(r — iacosd). 
For the azimuthal modes m = ±2 we get these expressions. 



(9) 



V27f 4 
And for the azimuthal mode m = we get 



o (r - 2M)(-3r + 8A/ + 3(r - 7M)Jl - ^) 
i{cose±lfLS 



(r(l + ./l-^)-M) 



2tt 



(r-2M)(-l + 2./l-if ) 



(r(l + Jl-if)-Af) 



(10) 



(11) 



^ Q (r- 2M)(-3r + 8M + 3(r- 7Af)./l - ^) 
-zsin^ 6'L5 



n/2^ 2 



r3(r(l + .>/l-^)-Af) 



(12) 



To all these expressions, we also need to add metric contributions. They have the same form as the ones in 
therefore we will not list them here. 



C. Initial data for the Teukolsky function — Case II 



Again, following the construction as outlined in the last section, and using the close-slow approximation, we arrive 
at the following expressions for the initial spatial metric and extrinsic curvature: 



Kab 



SSL 
Rr^ 



2 cose 

r-2M 
sin 6(4: cos^ 1) 



sine(4cos ^ — 1) 



— sin 26 sin 6 sin 2cp 



y/l-2M/r 

- sin 20 sin 6* sin 2o5 r • 
■ ^ ^smt'fl 



■yi~2M/r ^l-2M/r 

2r COS cos^ iy9 I sin0(l + cos^ 6*) sin2i^ 
3s^ 6*) sin 2iy9 — I sin sin 26* sin 2(y9 



(13) 



The metric again has a form, identical to the Misner solution 
in ||], the initial data for the Teukolsky function, is: 
For the azimuthal modes m = ±2 we get. 



Using the methodology and expressions we discussed 



3, ,2 (^-2M)(-l + 2Jl-2M) 

27r 4 , /t 2M^ 



(14) 



1 - ^) - M) 



J2i 



Q (r - 2Af)(-3r + 8Af + 3(r - 7Af) Jl - ^) 
T-(cos6i± 1)2^5 



(15) 



1 „ Ml) „ M) 



And for the azimuthal mode, m = the initial data are identically zero. To all these expressions we need to add 
the metric contributions, that have a form identical to the ones in our past work p|. 
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III. EVOLUTION OF THE DATA USING THE TEUKOLSKY EQUATION 



Given the Cauchy data from the last section, the time evolution is obtained from the Teukolsky equation 



a'^ sm^ - : — 



AMar ^ 
Ott + 4 



r + la cos a — 



A 



A^a, (A-^a,) + -—de {smOde) + 



1 



sm 



a 
A 



a{r — M) i cos ( 



{Acot'^9 + 2) ^^- = 0, 



(16) 
(17) 



where M is the mass of the black hole, a its angular momentum per unit mass (which is zero in our case), S = 



+ a cos 9, and A 



2Mr + a . Evolving the initial data we just calculated, with this equation will enable us 



to extract gravity wave waveforms that correspond to the late stage merger of spinning holes. We can also estimate 
the energy carried away by these gravitational waves. The radiated energy is given by [[l2| , 



dE 
'dt 



lim 

r — ^oo 



dn 



t 2' 

di ^{t,r,9,if) 



, dD, — sin 9 dd dip. 



(18) 



Note that the angular momentum radiated in the cases we consider is zero. 



IV. RESULTS OF THE EVOLUTIONS 



In this section we show waveforms and plots for energy radiated from the collision of two spinning holes. Recall 
that the two holes have equal mass and equal and opposite spin, and we are considering two cases. Case I in which 
the spins are aligned along the common axis of the two holes, and Case II in which the two spins are perpendicular 
to the common axis of the two holes. 

The waveforms that follow, are for a collision of two black holes that were initially separated by a conformal distance 
of 0.91 and had an individual spin of 0.1 in units of ADM mass. 

In figure we show the m = 2 mode of the Teukolsky function for a Case / collision. We see the typical quasi- normal 
ringing, in both the real and imaginary parts of the function. In figure |^ we plot the m — 2 mode of the Teukolsky 
function for a Case //collision. Note that the initial data is purely real for this case, hence the imaginary part of the 
waveform is identically zero. Quasi-normal ringing is self-evident here as well. This is the type of signal that gravity 
wave observatories like LIGO, will detect if they happen to witness a collision of the kind we are considering. 

Let us turn now to the evaluation of the radiated energies for these two cases. Figure |^ shows the radiated energy 
as a function of the initial spin, for a fixed separation of the holes. The first noteworthy thing is that the two curves, 
corresponding to Case / and Case II are very close to each other throughout. This indicates that the geometry of the 
initial configuration does not matter much as far as energy loss is concerned. But, it should be observed that for any 
spin S, Case / radiates more than Case II. The next important observation (by looking at the trend of the curves for 
large spin) is that, even for high values of the initial spin of the two holes this type of collision radiates less than 1% 
of its total mass. This means that inclusion of spin does not dramatically change our original estimate ||ll[ of a 
percent of the mass-energy being carried away by gravity waves in binary black hole inspiral. 



V. CONCLUSIONS 



We performed a first order, perturbative calculation to study the merger of two spinning holes. The exhaustive 
cases considered had holes with spins either aligned or perpendicular to the common axis of the two holes. For these 
cases we extracted waveforms and estimated energy radiated. We observed that in the case in which the spins of the 
two holes are perpendicular to the common axis, the system radiates a little less than the case in which the spins are 
along the common axis. In either case, the energy radiated appears to be a fraction of a percent of the total mass of 
the system. Since we know that the late stage inspiral of two non-spinning holes typically radiates about a percent of 
the mass-energy [|l^ Q this implies that adding spin to the holes doesn't really change such estimates by much. 

With these results (along with our earlier ones ||^) we can predict the properties of the gravity waves radiated from 
the ringdown phase of the slow inspiral of two spinning, equal mass black holes. 
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FIG. 1. This figure depicts the m = 2 mode of the Teukolsky function for a Case I collision. The solid line depicts the real 
part whereas the dotted one, the imaginary part. All quantities are in units of ADM mass. 
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FIG. 2. This figure shows the m = 2 mode of the Teukolsky function for a Case II collision. The solid line depicts the real 
part whereas the dotted one, the imaginary part. All quantities are in units of ADM mass. 



5 




FIG. 3. The radiated energy in a head-on collision of two spinning black holes as a function of the individual initial spin 
(that is equal and opposite for the two holes), for a fixed conformal separation of 0.91. All quantities are in ADM mass units. 
Looking at the behavior of the curves for high spin, one can estimate that less than 1% of the mass of the system is lost in a 
typical collision involving spinning holes. 
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